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ABSTRACT 

We examined the physical origin of the colour-magnitude relation of elliptical galaxies 
in clusters of galaxies in the hierarchical clustering scenario by using a semi-analytic 
model of galaxy formation. The semi-analytic model includes some physical processes 
connected with the galaxy formation, such as merging histories of dark haloes, gas 
cooling, star formation, supernova feedback, mergers of galaxies, and so on. There- 
fore it is easier to investigate how the galaxy formation and evolution depend on 
these physical processes by using the semi-analytic model rather than other numerical 
methods such as N-body and hydrodynamical simulations. Since the galaxy formation 
process is very complicated, analysis by using the semi-analytic model is effective. 
We particularly investigate the effects of the star formation process, the supernova 
feedback and the UV background on the formation and evolution of elliptical galaxies. 

It is suggested by observation that the colour-magnitude relation reflects the 
metallicity-luminosity relation of galaxies, and that for reproducing the colour- 
magnitude relation it is important to suppress the chemical enrichment of stars in 
faint galaxies. For suppressing the chemical enrichment, we propose some solutions. 
One is assuming a strong supernova feedback. Another is introducing the effect of 
the UV background radiation. These solutions degenerate, but they may be distin- 
guished by studying other properties of galaxies, like the luminosity function, colour 
distribution, etc. 

Moreover, we find that the supernova feedback model in the starburst which make 
the bulge component is important when the star formation time-scale in disc is longer 
than the age of the Universe. On the other hand, when the star formation time-scale 
is short, the resulting colour-magnitude relation does not depend on such a starburst 
model because most of stars are formed not by the burst but by the ordinary star 
formation process in disc. 
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1 INTRODUCTION 

It is well known that elliptical galaxies in clusters of galax- 
ies have the tight correlation between their colours and 
magnitudes (so-called 'colour-magnitude relation'; hereafter 
CMR). For example, the rms scatter about the mean CMR 
is typically ~ 0.04 mag in Virgo and Coma clusters of galax- 
ies. It is a comparable size to observational errors (Bower, 
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Lucey & Ellis 1992). Because it is believed that this tight 
relation reflects the formation and evolution processes of el- 
liptical galaxies, many people have studied this relationship 
in order to understand the galaxy formation process. 

The traditional scenario of the formation of ellipticals 
is that a monolithic protogalactic cloud collapses and then 
forms stars for a short time-scale until blowing the galactic 
wind (Larson 1974; Arimoto & Yoshii 1986, 1987). In this 
framework, the conventional interpretation of the CMR of 
ellipticals is considered as follows. Colour of each galaxy re- 
flects its mean stellar age and metallicity, while luminosity 
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of each galaxy reflects its mass. Therefore, there are two 
scenarios to make the difference of colours as follows. One 
is that redder luminous galaxies are formed at earlier stages 
of cosmological structure formation, and bluer faint galaxies 
are formed recently. Another scenario is that elapsed time 
of redder galaxies between the beginning of star formation 
and the epoch of galactic wind is longer than that of bluer 
faint galaxies, then more metal-rich stars are formed in lumi- 
nous galaxies. Therefore the luminous galaxies have redder 
colour. These two effects, age and metallicity, for reddening 
the luminous galaxies degenerate (so-called 'age-metallicity 
degeneracy'; Worthey 1996 and Arimoto 1996). This degen- 
eracy is solved in the collapse/wind model by Kodama & 
Arimoto (1997) by comparing the CMRs of high redshift 
cluster galaxies with those of theoretical models. If the CMR 
reflects the age sequence, young faint galaxies should be very 
blue at high redshift, and the slope of the CMRs evolves 
rapidly. However, the difference between the slopes of the 
CMRs at the present epoch and those at high redshift is 
small compared to the theoretical prediction in the age se- 
quence. This small difference is in agreement with the pre- 
diction in the metallicity sequence. Moreover, they suggested 
that age of ellipticals is longer than ~ 10 Gyr. Thus, they 
found that the CMR reflects the metallicity sequence. 

In contrast to such a collapse/wind model, recent devel- 
opments on both theory and observation of the cosmological 
structure formation are revealing that objects such as galax- 
ies and clusters of galaxies are formed through hierarchical 
clustering of smaller objects. Kauffmann & Chariot (1998; 
hereafter KC) applied their semi-analytic model to the prob- 
lem of the CMR, and reproduced the observational proper- 
ties of the CMRs, when the model includes the chemical 
evolution process and strong feedback to interstellar media 
by supernovae. Moreover, they found a tight luminosity- 
metallicity relation. 

KC interpreted why the CMR is reproduced when the 
feedback process is strong as follows. They considered that 
bright ellipticals must be formed by mergers of larger spiral 
galaxies but not from progenitors of faint ellipticals. Be- 
cause the feedback strength depends on the mass of galax- 
ies, the metallicity of stars in more massive progenitors be- 
comes higher than that in small progenitors. Therefore mas- 
sive ellipticals become redder than small ellipticals. How- 
ever, in the hierarchical clustering scenario, larger objects 
are formed by mergers of smaller objects. So such massive 
progenitors are formed by mergers of smaller galaxies. Thus 
the above conclusion may be insufficient for accounting to 
make the CMR, and it is needed to investigate how the 
metal abundance evolves in galaxies. The aim of this pa- 
per is, therefore, exploring what process, such as feedback 
and star formation, affects the CMR and the evolution of 
metallicity. Especially we focus on the slope of the CMR. 

In Section 2, we describe the semi-analytic model used 
here brieffy. In Section 3, we show the colour-magnitude re- 
lation in almost the same model as that of KC, and show 
evolutions of mean stellar metallicities of galaxies. In Sec- 
tion 4, we investigate the parameter dependence and some 
effects on the CMR. Section 5 is devoted to conclusions and 
discussion. 



2 MODEL 

In this section, the semi-analytic model which we use is de- 
scribed. Since the semi-analytic model includes some phys- 
ical processes connected with the galaxy formation, it is 
easy to understand the galaxy formation process, in con- 
trast to N-body and hydrodynamical simulations. Therefore 
we must include the following ingredients in our models at 
least: merging history of dark haloes, cooling and heating 
processes of baryonic gas, star formation and feedback by 
supernovae, mergers of galaxies, and stellar population syn- 
thesis model. The following procedure is mainly based on 
Kauffmann, White & Guiderdoni (1993), Cole et al. (1994) 
and Somerville & Primack (1998). 

The outline of the procedure is as follows. At first, the 
merging paths of dark haloes are realized by the extension 
of the PS formalism. Next, in each merging path, evolu- 
tion of the baryonic component, namely, gas cooling, star 
formation, and supernova feedback, are calculated. We rec- 
ognize a system consisting of the stars and cooled gas as a 
galaxy. When two or more dark haloes merge together, there 
is a possibility that galaxies contained in progenitor haloes 
merge together. The judgment whether the galaxies merge 
together or not is based on the dynamical friction time-scale. 
Here we define the lifetime of haloes as the elapsed time be- 
tween the formation time of the halo and the time at which 
the halo is subsumed into a larger halo. If the dynamical fric- 
tion time-scale is shorter than the lifetime, the galaxies will 
merge. If not, the galaxies will not merge and the common 
dark halo has two or more galaxies. Finally, we calculate the 
colour and luminosity of each galaxy from the star forma- 
tion history of each galaxy. Through the above procedures, 
we obtain a colour-magnitude diagram of galaxies. 



2.1 Merging history of dark haloes 

The merging histories of dark haloes are realized by using 
the method developed by Somerville & Kolatt (1999). They 
improved the model by Kauffmann & White (1993), based 
on the extension of the Press-Schechter formalism (Press & 
Schechter 1974) by Bower (1991), Bond et al. (1991) and 
Lacey & Cole (1993), for satisfying the mass conservation of 
haloes. 

Their algorithm is simple and can be quickly performed 
as follows. 

First, we consider two density fluctuation fields with 
two different smoothing scales, which corresponds to mass 
scales Ml and Mo (Mi < Mo). The variances of the density 
fluctuation flelds are ^i and 5*0, respectively. The variance Si 
corresponds to mass scale Mi uniquely if we fix a power spec- 
trum of the density fluctuation fleld. The linear density con- 
trast collapsing just at Zi is obtained as Ui — 1.69(1 -I- Zi) by 
the spherically symmetric collapse approximation (Tomita 
1969; Gunn & Gott 1972). 

Next we consider the probability distribution function 
of AS = Si - So, p{AS,Au})dAS, whore Auj ^ coi - ujq- 
When the density fluctuation flelds are Gaussian random 
flelds, this probability is simply 



p{AS,Auj)dAS ^ 



1 



Auj 



2^ (AS)3/2 



exp 



(Ao^ 
2AS 



dAS. (1) 



This function is a simple Gaussian distribution, if we change 
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the variable to a; = Alo/V^S. After this, we interpret the 
variance S as the mass M. 

The following procedure is based on a Monte Carlo 
method. First we set a mass Mo, or 5*0, as a mass of a finally 
collapsing halo, e.g., a cluster scale or the 'Milky Way' scale, 
and 'time-step' Aoj. By using eq.([ll), we pick out AS ran- 
domly. From this AS, Si, that is. Mi, is determined. This 
Ml is interpreted as a mass of a progenitor halo of the halo 
with mass Mo, collapsing at z = zi estimated by Aui. If 
Ml < Ml , where Mi is a minimum mass for an object iden- 
tified as an isolated halo at a corresponding redshift, the 
mass Ml is added to a diffuse 'accretion mass'. 

The above process is repeated until the rest of mass of 
the halo becomes less than Mi. In the same way, merging 
paths are recursively realized by regarding the halo with 
mass Ml as a starting point halo such as the halo with mass 
Mo. 

Mass functions at various redshifts obtained by using 
the above procedure reproduce the Press-Schechter mass 
function well. Of course, the Press-Schechter formalism itself 
has some problems, e.g., this formula does not include the 
spatial correlation of density fluctuations (Yano, Nagashima 
& Gouda 1996; Nagashima & Gouda 1997). Porciani et al. 
(1998) consider the effect of the spatial correlation in their 
method for realizing merging paths approximately. However, 
in spite of many theoretical uncertainties, the difference be- 
tween the Press-Schechter mass function and the mass func- 
tions given by N-body simulations is only a few factor in 
the number of haloes. This curiosity is reviewed by Monaco 
(1998). Therefore, for simplicity, we use the extension of the 
PS formalism in this paper. 

In this paper, we realize merging paths of progenitors 
of haloes with mass 5 x 10^ Mq , which corresponds to the 
circular velocity 14 ~ W^ km s~^. This is a typical size of 
clusters of galaxies. The time-steps are fixed to Az — 0.05. 
The minimum mass for identifying as an isolated halo is 
Ml = 10^°M©. Smaller objects than Mi may not be able to 
collapse because the mass of such objects are smaller than 
the Jeans mass (~ lO^^M©) after the reionization of inter- 
galactic medium by the UV background radiation at 2 ~ 10 
(Ostriker & Gnedin 1996), and because galaxies are formed 
actively at 2 ^ 5. We find that results presented in this paper 
does not depend on the choice of Mi ^ IO'^'^Mq. The power 
spectrum P{k) of density fluctuations given by Bardeen et 
al. (1986) is adopted, 



P{k) (xkr{k). 



2.34g 



1 + 3.1 



' + (16.1g)^ + (5.46g)^ + (6.71g)*] 



-1/4 



k 



(2) 

(3) 
- (4) 

where Q is the cosmic density parameter and h is the Hub- 
ble parameter. Ho = lOOh km s~^ Mpc"'^. We adopt the 
standard CDM model, f7 = 1 and h = 0.5, in this paper. 

2.2 Gas cooling 

In each realized merging path, we evaluate the amount of 
cooled gas through the radiative cooling, that is a direct 
material for forming stars. 



When a dark halo under consideration collapses, we es- 
timate the mass of hot diffuse gas contained in the halo. 
The sources of the diffuse gas are hot gas contained in its 
progenitor haloes and in the accreting matter. The baryon 
fraction of the accreting matter is defined as Qb/i^, where Qb 
is the cosmic mean baryon density. The density distribution 
p(r) of this hot gas in the halo is assumed as the isothermal 
distribution, p(r) ex r~^. The absolute value of the den- 
sity profile is determined by mass and collapsing redshift of 
haloes under the spherical collapse model. 

The cooling time-scale TcooI is obtained as a function of 
the radius from the hot gas density profile, the temperature 
of the hot gas, and the cooling function A(r) as follows, 

where iinip is the mean molecular weight, n^ir) is the elec- 
tron number density at the radius r, and k is the Boltzmann 
constant. The cooling function depends on the metallicity of 
the hot gas. We adopt the both cooling functions of zero- 
metallicity and solar-metallicity given by Sutherland & Do- 
pita (1991), and determine the cooling efficiency depending 
on the metallicity of the gas from the two cooling functions 
by interpolation and extrapolation. 

When a halo merges with another halo of compara- 
ble size (when the mass ratio of these haloes is larger than 
f reheat), the diffusc hot gas contained in the halo is shock 
heated to the virial temperature of the halo. In this case, the 
hot gas within the cooling radius Vcooi cools and the cooled 
gas mass is simply added to the cold gas reservoir of the 
galaxy. The cooling radius is a radius at which the cooling 
time-scale equals to the time-step At of the merging histo- 
ries of dark haloes. On the other hand, when the halo merges 
with small haloes whose mass ratio to the halo is less than 
f reheat, the hot gas is uot shock heated and preserves its 
temperature, and the density profile is not disturbed. This 
parameter is introduced by Somerville & Primack (1998) 
firstly. In such a case, gas contained in the smaller haloes 
may remain in outer region of the new halo. Therefore, we 
introduce the following process when the mass ratio of the 
merging haloes is less than f reheat- Now we define teiapse as 
the elapsed time since the last shock heated epoch of the 
halo and rcooiifeiapse) as a radius at which the cooling time- 
scale equals to the elapsed time. In this case, because the 
density profile is maintained and the gas within reooi{teiapse) 
has already cooled, the cooled gas mass is estimated as 

AMeooi = / /9(r)47rr^dr. (6) 

This procedure is adopted only when reooi{teiapse)/R < 0.8. 
When rcooi/R > 0.8, the hot gas may fall into the centre of 
haloes and cool because most of hot gas has already cooled 
at the previous time-step. So the gas density profile is recal- 
culated and rcooi is estimated by TcooI = tufe as in the case 
of the major merger of haloes, while the density profile of 
the dark halo will be preserved. 

We can avoid the dependence of the amount of the 
cooled gas mass on the artificial time-step of calculation by 
using the above procedure. In this paper, we use f reheat ~ 
0.2. Of course, if we use the block model which was used by 
the Durham group (Cole & Kaiser 1988; Cole et al. 1994), 
this parameter is not required. In the extension of the PS 
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formalism, masses of hafoes are given by randomly and the 
time-step is fixed. On the other hand, in the block model, 
density contrasts which determine the time-step are given 
randomly, while masses of haloes are divided by a factor of 
2. So mass of a new collapsing halo is always twice of a max- 
imum progenitor. Thus the mass ratio of the merging haloes 
are usually larger than freheat = 0.2. 

In the case that the circular velocity of haloes exceeds 
500 km s~^, we prevent the cooling process by hand to sup- 
press the formation of extremely big monster galaxies. If we 
consider only the mean density of the hot gas, the maxi- 
mum mass of galaxies is ~ IO^^Mq (Rees & Ostriker 1977). 
In our calculation, since we adopt the isothermal distribu- 
tion p{r) oc r~^ , the cooling process becomes efficient in 
the central region of haloes. Actually this phenomenon is 
suggested observationally (e.g., Fabian 1994). However the 
number density of bright galaxies is much larger than the 
observed one, unless we prevent the cooling process in such 
large haloes. Thus we adopt the same unsatisfactory solution 
as Kauffmann et al. (1993). Somerville & Primack (1998) 
also adopt the same solution but for 400 km s~^ of the cut- 
off circular velocity. A detailed study is needed in order to 
understand this effect physically. 



2.3 Star formation 

In Kauffmann et al. (1993), the star formation in disc is 
described by the following simple law, 

M. = ^^^, (7) 

r. = r°(l + 2)-'/^ (8) 

where eq.(0) denotes the rate of stars newly formed. This 
form is the same as the one adopted by Kauffmann et al. 
(1993), and referred as the 'Munich model' by Somerville & 
Primack (1998). The redshift dependence on the star forma- 
tion time-scale r« is obtained as follows. We adopt the simple 
star formation interpretation, r, ex R/Vc, that is, r* is pro- 
portional to the dynamical time-scale. From this, the virial 
theorem and the spherical collapse model, we obtain that 
r* ex (1 -f z)~^'^. The value of the star formation time-scale 
at 2; = 0, r°, is a free parameter. This parameter describes 
the star formation with a long time-scale in a disc of spiral 
galaxies. Moreover we need to fit the amount of the cold 
gas in the 'Milky Way'-like spiral galaxies, which have the 
same circular velocity as our Galaxy (Milky Way), Vc — 220 
km s~^, to the amount of the cold gas in our Galaxy. If r° 
is small, almost all of the cold gas is transformed into the 
stars until the present epoch. So rf is larger than a few Gyr. 
Stars formed by the above process constitute a disc of 
galaxies. Star formation in a bulge occurs when two or more 
gala xies with comparable masses merge together (see Sec- 
tion 2.5). In this merging process, since the system is dy- 



namically disturbed, starburst may occur. Then both stars 
formed in the burst and disc stars may fall into the central 
region of the galaxy and constitute the bulge. We describe 
this event by adopting very short star formation time-scale. 
Of course, the above star formation laws are oversim- 
plified. In the present situation, stars are considered to be 
formed as follows. The star formation process may be de- 
termined by the local properties, such as density and tem- 



perature. If we consider only atomic cooling in our calcula- 
tion, the gas cools until ~ 10* K. So we identify such gas 
cooled until ~ 10* K as the cold gas. However, the pres- 
sure of the gas is still too high to form stars. The cooling 
by molecules is needed to cool the gas lower than 10* K. 
Then the gas cools until a few K by the molecular cooling, 
and fragments into small clumps. These clumps evolve to 
stars. These processes are performed in molecular clouds, 
whose mass is about IO^Mq. Therefore we should consider 
the processes of the molecular cooling and the fragmentation 
of the 'cold' gas (~ 10* K). It is studied by considering the 
ensemble of molecular clouds to connect the star formation 
in molecular clouds with the global star formation properties 
in galaxies (Fujita 1998; Fujita & Nagashima 1999). These 
processes are very complicated, so we do not adopt them 
in our semi-analytic model, since the simplicity which is an 
important merit of this model will be lost. 



2.4 Supernova feedback 

When stars form, massive stars with short lifetime explode 
and heat up the surrounding cold gas. The number of super- 
novae per solar mass of stars formed is r]sN — 7 x 10~'^Mq^ 
for the Salpeter IMF (Salpeter 1955). The kinetic energy 
from a supernova, Esn, is about 10^^ erg. When a fraction 
e of the released energy is used to heat up surrounding cold 
gas, the amount of this reheated gas is evaluated as 



l\Mrehe.at — £- T72 '^*- 



(9) 



This feedback process has many uncertainties actually, and 
therefore we adopt a simple description (Cole et al. 1994), 



AMr, 



\Vhot) 



ALAt = i3PLAt, 



(10) 



where Vhot and ahot are free parameters. Eq.(y), which is 
adopted by Kauffman et al. (1993) corresponds to Uhot ~ 2, 
and Cole et al. (1994) used Uhot = 5.5 in their fiducial model. 
It is the conclusion of KG that this feedback process is a 
key parameter to determine the slope of the GMR. Therefore 
we also investigate the dependence of GMR on the feedback 
process. 

2.5 Mergers of galaxies 

When two or more haloes merge together, the new common 
halo may have two or more galaxies. In that case, the galax- 
ies will lose their energy due to dynamical friction, and then 
fall into a centre of the new common halo. Finally they may 
merge together as the case may be. 

When a halo collapses, a central galaxy contained in the 
largest progenitor is identified as the central galaxy of the 



new common halo. The cooled gas described in Section 2.2 
accretes to this central galaxy. Other galaxies are identified 
as satellite galaxies. 

We calculate the elapsed time from the epoch at which 
each satellite galaxy was identified as satellite, namely, it was 
contained in a common halo not as a central galaxy but as 
a satellite galaxy. It should be noted that this elapsed time 
is reset to when the 'major merger' of dark haloes occurs 



(see Section 2.2) because orbits of satellite galaxies may be 
violently disturbed by the major merger of haloes. When 
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this elapsed time exceeds the dynamical friction time-scale 
for the galaxies in the new common halo, the satellite merges 
with the central galaxy. The dynamical friction time-scale is 
(Binney & Tremaine 1987) 



Star 



"""^^ - InAGM,, 



(11) 



where R and Vc are the virial radius and circular velocity of 
the new common halo, respectively, M^at is the total mass 
of a halo to which the satellite belonged as central galaxy, 
and In A is the Coulomb logarithm, which is approximated 
as ~ ln{l + Mff / M^^i) (Mh is the mass of the common halo) 
(Somerville & Primack 1998). 

In this paper, we do not treat mergers between satel- 
lite galaxies for simplicity, which is described in Somerville 
& Primack (1998). There is a possibility that this effect 
changes the resulting CMRs. This effect will be considered 
in future work. 

When a satellite galaxy merges with a central galaxy, 
and moreover the mass ratio of the galaxies is larger than 
f bulge, all stars of the satellite galaxy and disc stars of the 
central galaxy are incorporated with the bulge of the cen- 
tral galaxy. Then cold gases of both galaxies turn to stars in 
the bulge with a very short time-scale ('starburst'). In this 
paper, we use the following two burst model on the feed- 
back process. One is that the strength of the feedback is the 
same as that in the disc, according to eq.([lO|) (model sbA). 
Another is that there is no feedback effect during the burst, 
that is, all cold gas turns to stars through the 'closed-box' 
process (model sbB). In order to decide which model is cor- 
rect, detailed simulations of merging process will be needed. 
On the other hand, if the mass ratio of the galaxies is smaller 
than f bulge, all stars of the satellite galaxy are incorporated 
with the disc of the central galaxy and the cold gas reser- 
voir of the satellite merges with that of the central galaxy, 
and then disc stars are formed from the cold gas with the 
star formation time-scale obtained from eq.(@). In this pa- 
per, we adopt ftuige = 0.2. This value is the same as that in 
Somerville & Primack (1998). They adopt this value for re- 
producing observational fractions of the number of galaxies 
of each morphology. 

2.6 Chemical evolution 

Chemical evolution is treated in almost the same way as 
described in KC. The instantaneous recycling approxima- 
tion is adopted (Tinsley 1980). We show the gas and metal 
exchange in Fig.tl schematically. Metals ejected from super- 
novae are recycled as follows. 

When two or more haloes merge together, metals con- 
tained in the hot gas is joined together. If two or more galax- 
ies merge together, metals contained in the cold gas reser- 
voirs of galaxies are joined together. When the hot gas cools, 
metals contained in the cooled gas are also joined to the 
cold component. Here we assume that metals in cold and 
hot components are well mixed in each component. Stars 
formed from the cold gas have initially the same metallicity 
as the cold gas. This metallicity affects colours of galaxies 
directly. 

The amount of metals ejected from supernovae is char- 
acterized by y, which is heavy element yield for each gener- 
ation of stars. The fraction / of the ejecta is ejected directly 




formation 



hot gas 



radiative cooling 



cold gas 



Figure 1. Schematic description of gas and metal recycling in 
star formation. / is a fraction of ejected mass from evolved stars 
into hot gas directly. 



into the hot gas, and the rest, 1 — /, is incorporated with 
the cold gas. We adopt / = 0.3, which is the same value as 
in KC. The gas fraction returned by evolved stars R is 0.75 
in this paper. Simultaneously, the supernovae heat up the 
surrounding cold gas, therefore metals contained in the cold 
gas are also returned to the hot gas. 



2.7 Stellar population synthesis and identification 
of morphologies 

In order to compare our results with observations such as 
CMRs directly, stellar population synthesis model must be 
considered. We adopt a model by Kodama & Arimoto (1997) 
as such a population synthesis model. Luminosity and colour 
of stars cannot be theoretically estimated without uncertain- 
ties because the stellar evolution models themselves include 
the uncertainties. For example, as for the convection of stel- 
lar gas, we have only a phenomenological theory ('mixing 
length theory'). However, their model is adequate for our 
purpose of investigating the slope of the CMR. Once we un- 
derstand the properties of galaxies such as the CMR physi- 
cally and qualitatively, even if the stellar population model 
is changed, we fit immediately our results with observations 
by changing parameters mentioned above. 

The IMF which we adopt is the Salpeter type with a 
slope of 1.35, and the mass range is O.IM0 ~ 6OM0. The 
range of stellar metallicity Z, of simple stellar populations 
is 0.0001 ~ 0.05. 



In Section 2.3, we divide the stellar component into disc 



and bulge components. Morphology of each galaxy is deter- 
mined by the B-band bulge-to-disc luminosity ratio (B/D). 
Simien & de Vaucouleurs (1986) showed that the Hubble 
type of galaxies correlates with the B-band luminosity B/D. 
In this paper, galaxies with B/D > 1.52 are identified as el- 
lipticals, 0.68 < B/D < 1.52 as SOs, and B/D < 0.68 as spi- 
rals, according to their results. It is shown that this method 
for classification reproduces observations well by Kauffmann 
et al. (1993) and Baugh, Cole & Frenk (1996). 



3 ORIGIN OF COLOUR-MAGNITUDE 
RELATION 

In this section, we explore the origin of the CMR. In this 
paper, because we investigate how the CMR depends on the 
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Table 1. 


Parameters for models. 








Model 


Vhot 


(km s 


-^) 


C^hot 


r° (Gyr) 


UV 


burst 


A 




280 




2 


20 


off 


sbA 


B 




100 




2 


20 


off 


sbA 


C 




280 




2 


20 


off 


sbB 


D 




100 




2 


20 


off 


sbB 


E 




280 




2 


2 


off 


sbB 


F 




100 




2 


2 


off 


sbB 


G 




240 




5.5 


2 


off 


sbB 


H 




100 




5.5 


2 


off 


sbB 


I 




400 




2 


2 


off 


sbB 


J 




280 




2 


2 


on 


sbB 



physical processes such as star formation, supernova feed- 
back, and so on, we fix the cosmological parameters to the 
standard CDM model, that is, n = 1, A = 0, // = 50km s"^ 
Mpc~^, Q,h ~ 0.06, and erg = 0.67. We refer to the models 
considered by the models A, B, C and D as shown in Table 
1. In this table, we also show models from E to J, which 
are adopted in the next section. The fifth column, UV, is 
explained in Section 4.3. The sixth column, burst, is shown 
in Section E^. The yield y is equal to 0.038 — 2Zq in all 
the models. In KC, y = 1.2Zq in the low feedback model. 
However, in order to see the effects of the feedback and so 
on, we fix the value of the yield in all models. 

In the following sections, we take notice of only the 
slopes of CMRs, because the luminosity of galaxies can be 
translated by considering the following reason. Stars are 
formed according to the IMF. Mass of luminous stars is 
larger than ~ O.O8M0, which is determined by the crite- 
rion of nuclear burning. However, there is a possibility that 
invisible stars with mass smaller than O.OSM© are formed. 
If there are many invisible stars in galaxies, the galaxies 
become faint compared to the case that all stars are lumi- 
nous. The ratio of the invisible stars to the luminous stars is 
treated as a free parameter in the previous work. Therefore 
the absolute value of luminosity can be adjusted by this pa- 
rameter. The suitable value of this parameter will be deter- 
mined by considering other observational quantities. Thus 
this parameter does not affect the slope of the CMR. 



3.1 Colour— magnitude relations 

In Fig.Q, we show the CMRs in the models A and B with 
Uhot ~ 2 which is the value adopted by KC. The dots denote 
galaxies identified as ellipticals, and the solid lines show the 
observational CMRs (Bower et al. 1992). The criterion to 
pick out ellipticals among all galaxies is shown in Section 
2^ 

As in KC, in the model A with high feedback efficiency, 
the slopes of the CMRs are in roughly agreement with the 
observations, but the dispersion is larger than that of obser- 
vations, ~ 0.04 mag. This is due to the recent star formation 
owing to long star formation time-scale, r^ = 20 Gyr. On the 
other hand, in the model B with low feedback efficiency, the 
slopes of the CMRs are nearly flat. At the bright-end, the 
colour becomes redder about 0.1-0.2 mag, and at the faint- 
end {Mv ~ —18), about 0.2-0.3 mag, compared to that in 
the model A. In order to see the physical relation between 
the reddening and the feedback intensity, we investigate the 
age- and metallicity-luminosity relations. 
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Figure 2. Colour-magnitude diagram of model A and B. (a) My 
v.s. U — V in the model A. (b) Same as (a), but for the model 
B. (c) My v.s. U — V in the model A. (d) Same as (c), but for 
the model B. The models are characterized in Table 1. The dots 
denote galaxies identified as ellipticals. The solid lines show the 
observed CMRs (Bower et al. 1992). 



We define the V-luminosity weighted mean stellar age 
and metallicity as follows, 






Jg° Lv{t, Z,ald)M,{t)Zcold{t)dt 

f^° Lv{t,Z,M)M,{t)dt 
Jl" Lv{t,Z,,M)M,{t)tdt 
J^"Lv{t,Zcoid)M*{t)dt' 



(12) 
(13) 



where Zcoid is the metallicity of cold gas, M*{t) is the star 
formation rate at t, and Lv{t,Zcoid) is the luminosity of 
stars at present, t = to, with age to — t and metallicity 
Zcoid- In both models, the distributions of mean stellar age 
are alike (upper panels of Fig.H), and the mean stellar ages 
spread about 5 ~ 11 Gyr. On the other hand, there are 
tight relations between the V-luminosity and the mean stel- 
lar metallicity. The slopes of the relations in the models A 
and B are almost the same, but the value of the metallicity 
in the model B is higher than that in the model A. In the 
model B, because ages of most of dwarf ellipticals are old 
and because the metallicity is high, the CMRs are fiatter 
than those in the model A. 

Next we show the CMRs in the models C and D, in 
which the all cold gas is transformed to stars in the star- 
burst process (FigJjt). In both models, the slopes of CMRs 
are fiat, as well as in the model B, while the number of galax- 
ies in the model C is less than that in the model D, owing to 
the high feedback efficiency in the model C. The difference 
between the models A and C is only in the feedback model. 
This difference will cau se t he difference in the metallicity 
of galaxies (see Section 3.3). In Fig.Bl We show the same 
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Figure 3. Upper panels: age-magnitude diagram, (a) Model A. 
(b) Model B. Lower panels: metallicity-magnitude diagram, (c) 
Model A. (d) Model B. The dots denote elliptical galaxies. 

figure as FigJS but for the model C and D. The age distri- 
butions of the models C and D are almost the same as those 
of the models A and B. On the other hand, the metallicity 
distributions are different. The mean metallicities of dwarf 
ellipticals increase by {Z,)ly ~ 0.03 and the distributions 
are flat in both models C and D, while {Zt)Lv ~ O-O^ ^^ 
Mv ~ — 19 in the model A. Therefore the difference of the 
feedback models mainly affects the metallicity of dwarf ellip- 
ticals. Thus we can say that the feedback model is essential 
to determin e the properties of the CMR, when r" — 20 Gyr 
(see Section 3.3). 

From these figures, we conclude that the CMR results 
from the metallicity-luminosity relation as KC concluded. 
The conclusion of the origin of the metallicity-luminosity 
relation in KC is that ellipticals are formed by mergers 
of gas-rich spirals, and that massive ellipticals are formed 
from massive spirals because the feedback efficiency of such 
massive galaxies is small, and so more metal-rich stars are 
formed (Fig. 2 in KC). We will show the evolutions of mean 
stellar metallicity and age of each galaxy in order to inves- 
tigate this process in detail. 



3.2 Evolution of mean stellar metallicity 

In this subsection, we show the evolutions of the mean stellar 
metallicity {Z*)a/, , which are weighted by stellar mass for 
simplicity. 

In the upper panels of Fig.H, we show the evolutions 
of stellar metallicity of luminous ellipticals picked out in 
the models A and B. Their final magnitudes and colours 
are My = -22.0, U - V = 1.47, V - K = 3.09 and 
Mv = -22.0, U-V = 1.61, V-K = 3.23, respectively. The 
open circles denote the mean stellar metallicities of progen- 
itors in the elliptical state, which is defined as the bulge- 
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Figure 4. Same as Fig.y, but for model C and D. Both models 
have a long star formation time-scale, t2 = 20 Gyr, and the 
starburst model sbB is adopted. Model C: V^ot = 280 km s~^. 
Model D: Vhot = 100 km s'^. 




Figure 5. Same as FigJsl but for model C and D. 

to-disc stellar mass ratio Altuige/Mdisc > 0.6, and the dots 
show those in the spiral state, Mi,uige/Mdisc < 0.6. These 
two luminous galaxies are identified as galaxies experiencing 
the same merging path of dark haloes. Because the merging 
histories of dark haloes of these galaxies are the same, the 
difference between results in different models is caused from 
only parameters shown in Table hi. These galaxies rapidly 
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grow in mass and the feedback effect becomes weak at early 
times. So the mean stellar metallicities also exceeds Zq at 
early times. Therefore, the difference between these galaxies 
is small. The final metallicity of elliptical galaxy in the model 
A (Fig.Ma) is smaller than the yield y = 0.038. The reason 
for this is considered as follows. One is that the feedback 
effect remains a little because of the shallower dependence 
of feedback on circular velocity, ahot = 2. Second is that the 
star formation time-scale is very long, r^ — 20 Gyr. On the 
other hand, the difference between faint ellipticals is large 
(lower panels of Figj6|). The ratio of the final metallicities 
is more than a factor of two. These faint ellipticals are also 
picked out in the same merging path. Their final magnitudes 
and colours are My = -18.4, U -V = 1.24, V - K = 2.86 
and Mv = -19.1,t/- V = 1.31,V-K = 3.07, respectively. 
In the high feedback efficiency model A, the growth of the 
mean stellar metallicity is slow and the final metallicity is 
very small, {^.)m, ~ 0.013. On the other hand, the growth 
in the low feedback efficiency model B is rapid and then the 
final metallicity is high, {Zf)M, ~ 0.028. This shows that 
in the low feedback model even small galaxies can evolve 
chemically. 

In this figure, some progenitors have periods of decrease 
of their metallicities. This is because gas with low metallicity 
accretes a galaxy with high mean metallicity and stars are 
newly formed from the gas. Therefore such stars have low 
metallicity and the mean metallicity of the galaxy decreases. 

In FigM we show the same figures but for the mod- 
els C and D with the starburst model sbB. The final 
magnitudes and colours of picked out luminous ellipticals 
are Mv = -22.2, t/ - V = 1.55,V ~ K = 3.19 and 
Mv = -22.0, U -V = 1.62, V - K = 3.25, in the mod- 
els C and D, respectively. Those of faint ellipticals are 
Mv = -18.8, [/ - V = 1.45, y - K ^ 3.13 and Mv = 
-19.2, U -V ^ 1.36, V -K ^ 3.14, respectively. In aU four 
panels, the evolutions of metallicities of galaxies in the spiral 
state are the same as those in the model A and B because 
the difference between the models A and B and the mod- 
els C and D is only the starburst model. In contrast, the 
metallicities of galaxies in the elliptical state are different, 
especially in the high feedback model and in the dwarf ellip- 
ticals (upper- and lower-left panels and lower-right panel in 
Fig.0) . Moreover, the final metallicities of these four galaxies 
are almost the same, {Z,)m, ~ 0.043. This value is nearly 
equal to the value of the yield y. This shows that the star- 
burst increases the metallicity to the yield, when all cold 
gas turns to stars with very short time-scale. We will ex- 
plained this increase of metallicity in the starburst in the 
next Section. 

From the above results, we conclude that the colour 
of luminous galaxies reflects the yield strongly because the 
effect of the feedback on such luminous galaxies is very 
weaker than that on faint galaxies. On the other hand, in 
faint ellipticals, their colour reflects the feedback efficiency, j3 
[— (F/Vhot)~"''°* , see eq.(px|)], directly, when the feedback is 
effective even in the starburst process. In other words, when 
the value of y is fixed, the evolution of stellar metallicity 
in faint galaxies is suppressed in the case of high feedback 
models, and the metallicity of luminous galaxies is almost 
determined by the yield, y, because the feedback effect be- 
comes negligible in such luminous galaxies. Therefore, a sup- 
pression of growth of stellar metallicity is required in small 
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Figure 6. Evolution of mass-weighted mean stellar metallicity. 
Each symbol shows the mean metallicity of a progenitor at the 
time. The open circles denote progenitors in the elliptical state, 
and the dots the spiral state. These states are defined in the text, 
(a) A luminous elliptical in model A. Mv = —22.0, U — V = 
1.47, V - K = 3.09. (b) A luminous elliptical in model B. My = 
-22.0, U-V = 1.61, V-K = 3.23. (c) A faint elliptical in model 
A. Mv = -18.4, U-V = 1.24, V-K = 2.86. (d) A faint elliptical 
in model B. Mv = -19.1, U -V = 1.31, V - K = 3.07. 



galaxies for making the slope of CMRs. The relationship 
between the metallicity and the feedback is discussed in the 
next section. 



3.3 Roles of feedback and yield 

In a chemical evolution of a simple monolithic cloud collapse 
model, the evolution of metallicity of cold gas for a constant 
star formation time-scale, r», is described as 



^cool{t} 



(1 - ^cooi)exp 



-(i-/)«y 



t-ts 



(14) 



where a is a locked-up mass fraction, a = 1 — _R (_R is the 
gas fraction returned by evolved stars) , and Z^^^^i is an initial 
metallicity of cold gas at i = t^. The rate of increase of stellar 
mass is 



M,{t) = (l-i?) 



M^oolit) 



(15) 



where we use the instantaneous recycling approximation. 
From these two equations, the mass-weighted mean stellar 
metallicity is obtained as 



{Z^{t)), 



^ ^ l-exp[-(l + /3-(fl-at/)(l-/))^] 

l-exp[-(l + /3-i?(l-/))i^] '^ 



where 
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Figure 7. Same as Fig.1 
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j, but for model C and D. (a) My 

- K = 3.19 in model C. (b) My 

- K = 3.25 in model D. (c) My 

- K = 3.13 in model C. (d) My 
K = 3.14 in model D. 
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Figure 8. Same as Fig. El but for model E and F. Both models 
have a short star formation time-scale, t2 = 2 Gyr, and the star- 
burst model sbB is adopted. Model E: Vhot = 280 km s~l. Model 
F: Vhot = 100 km s-i. 



F = 



^cool 



1 + 



(l-/)Qi; 



. (17) 

l+l3-R(l-f) 

When the star formation time-scale is enough smaller than 
the Hubble time or when the starburst occurs, the final 
metallicity is determined by F {tJT, — > oo). The form of F 
shows that when the feedback is strong the final mean stel- 
lar metallicity becomes small, {Zti^jr, -^ oo))a/, —>■ Z'^ooi 
as l3 —> oo. Moreover, increasing y is corresponding to de- 
creasing the feedback strength /3. This effect determines the 
absolute value of final mean stellar metallicity, especially for 
faint galaxies. On the other hand, when the feedback is neg- 
ligible (/3 ^ 1), the factor F is determined by only y. So the 
efTect of y is remarkable in luminous galaxies rather than 
in faint galaxies when the starburst model sbA is adopted, 
which is the model that the supernova feedback strength re- 
mains even in the starburst. On the other hand, when the 
starburst model sbB is adopted, which is the model that all 
the cold gas turns to stars with very short time-scale, the 
final metallicity of elUpUcals does not depend on their mass, 
because /3 = / = in the starburst. In this case, we find that 
the final metallicity (Z.)m. ~ 0.04 when Z°^„i ~ - 0.01. 
This value is almost the same as the final metallicities of 
ellipticals in Fig.n. Thus the results in the previous section 
are mostly explained by eqs.dlq) and ([L7|). 



4 PARAMETER DEPENDENCE 

In this section, we investigate the parameter dependence of 
the CMR on the star formation time-scale and the circular 



velocity dependence of the feedback. In Section 4.3, we also 
investigate the effect of the UV background radiation. 



4.1 Star formation time-scale 

In the previous section, the star formation time-scale, r° — 
20 Gyr, is longer than the Hubble time. Next we investigate 
how the CMR changes in the case of a short star formation 
time-scale. We adopt r^ = 2 Gyr in this section. 



4-1.1 CMR and starburst model 

We show the CMRs in the models E and F (Fig.|). The 
starburst model sbB is adopted in these models. 

In the left panels (the model E), the dispersion becomes 
small compared to the models A and C and the CMR repro- 
duces the observations well. The colours of ellipticals become 
redder compared to those of the model A (Fig.0) and those 
of bright ellipticals of the model C (Fig.l4|) , and become bluer 
than those of dwarf ellipticals of the model C. Because the 
starburst model of the model E is the same as that of the 
model C, these figures show that the star formation time- 
scale is one of the essential parameters for the slope of the 
CMR when we adopt the starburst model sbB. 

In the right panels (the model F), luminous ellipticals 
are very red and the CMR is flat at My ;$ — 20. The colour 
of the giant ellipticals nearly corresponds to the colour when 
Z, ~ 0.04 and the age of galaxies is nearly equal to 10 Gyr. 

Moreover, we find that the dependence of CMRs on the 
starburst model becomes negligible. The above results do 
not change when the starburst model sbA is adopted (not 
shown) . 

These properties are discussed in Section 4.1.2| . In the 
same way as Section pi we investigate the metallicity and 
age of ellipticals next. 
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4-1.2 Metallicity and age of ellipticals 

In order to understand the properties mentioned above, we 
show the metallicity- and the age-luminosity diagrams in 
Fig.H. Ages of ellipticals become older about 1-2 Gyr than 
those in the case of r^ = 20 Gyr on average. However, 
this difference corresponds to less than 0.1 mag in colour. 
We cannot explain the origin of the difference between the 
CMRs in the models from A to D and the models E and F 
only by the age difference. 

The metallicity in the model E [FigMc)] increases in all 
range of magnitude compared to the model A [Fig.g(c)], and 



{Z*)i 



0.028 at Mv 



-23 and {Z^)lv ~ 0.018 at Mv 



— 19. The reason of the metallicity increase is considered as 
follows. In eq.(hq) it is only the case of i/r, -^ oo that the 
final metallicity is determined by only F. This is correct 
when the star formation time-scale is quite shorter than the 
Hubble time. However, the value r^ = 20 Gyr is longer than 
the Hubble time and does not satisfy the above condition. 
The mass- weighted mean stellar metallicity {Z,{t))M, is a 
monotonousl y inc reasing function about t in the description 
in the Section 3.3. So the actual final mean stellar metallicity 
is lower than 1 ~ F. Therefore when we adopt a short star 
formation time-scale, mean stellar metallicities of galaxies 
increase totally and colours become redder, compared to the 
model A. 

In contrast to the model A, in the model C with the 
same starburst model sbB as the model E, the metallicity is 
higher than that in the model E. The difference between the 
metallicities in the models C and E is emphasized at dwarf 
elliptical scale. The reason that the metallicity of dwarf el- 
lipticals is low in the model E is considered as follows. In 
the spiral state, the cold gas turns to stars in very short 
time-scale, r, = 2(1-1- z)~^'^ Gyr. In this star formation 
mode, supernova feedback can reheat the cold gas, and most 
stars are formed in the spiral state owing to the short star 
formation time-scale. Therefore, when the major merger oc- 
curs, the number of stars formed through the starburst is 
small. Note that when stars are formed with the feedback, 
the stellar metallicity becomes lower than that in the case of 
the starburst without the feedback, /3 = / = (see Section 
3.3). Thus the metallicity becomes small compared to that 
in the model C, in which most of stars are formed during 
the starburst without the supernova feedback. Besides the 
difference between the feedback models during the starburst 
does not affect the metallicity because of the small number 
of stars formed during the starburst. 

Thus the metallicities of ellipticals in the model E are 
higher than those in the model A and lower than those in 
the model C, and the results do not depend on the feedback 
model. 



4.1.3 Evolution of mean metallicity 

Next we investigate the evolution of metallicities of some 
sample galaxies in order to show the time-scale of the metal- 
licity increase and that most of stars are formed in disc in 
the models E and F. 

We explicitly show the evolution of the mean metallicity 
of four ellipticals in FigJld. The magnitudes and colours of 
the picked out luminous ellipticals are Mv ~ —22.0, U — V = 
1.6l,V-K = 3.21 a.ndMv = -22.0,17-1/ = 1.77,V-K = 



(o) 
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Figure 9. Same as Fig.H, but for model E and F. 



3.34 in the models E and F, respectively. Those of the faint 
ellipticals are Mv = -18.5,(7 -V = 1.44, F - K = 3.07 
and Mv = -19.0, U ~V = 1.60, V -K = 3.24, respectively. 
This figure resembles FigJq of the models A and B with the 
starburst model sbA rather than Fig.M of the models C and 
D with the starburst model sbB. Note that the starburst 
model in the models E and F is the same as that of the 
models C and D. In the models A and B, there are many 
cases that the starburst occurs in the middle of the metallic- 
ity increase, while in the end of the metallicity increase in the 
models E and F. This shows that the star formation time- 
scale is short enough and that most of stars are formed in 
disc, where the supernova feedback is effective. In the mod- 
els C and D, many stars are formed during the starburst 
without the feedback, so that the mean stellar metallicity 
becomes very high (see Section 3.3). On the other hand, in 
the models A and B, the feedback is effective even during 
the starburst. So the resulting CMRs of the models E and 
F resemble those of the models A and B, while the mean 
metallicities of the models E and F are higher than that 
of the models A and B because of the short star formation 
time- scale. 

Thus most of stars are formed in disc and the differ- 
ence between the CMRs with the starburst models sbA and 
sbB is negligible when the star formation time-scale is short 
enough. 



4.2 Slope of colour— magnitude relation 

In this subsection we investigate how the CMR depends on 
the parameter ahot- We adopt ahot = 5.5 in the models G 
and H. This value of auot has been adopted by Durham 
group (e.g.. Cole et al. 1994). In the model G, Vhot = 240 
km s~^. In this case, the feedback strength (3 has the same 
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Figure 10. Same as Fig 
-22.0,(7 -V = 1.61, V 
-22.0,(7 -V = 1.77, y 
-18.5,(7 - y = 1.44, V 
-19.0,(7- y = 1.60, y- 



g but for model E and F. (a) My 

- K = 3.21 in model E. (b) My 

- K = 3.34 in model F. (c) My 

- K = 3.07 in model E. (d) My 
K = 3.24 in model F. 
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Figure 11. Same as FigJ3, but for model G and H. Both models 
have a short star formation time-scale, r^ = 2 Gyr, and a steep 
dependence of the feedback on the circular velocity. Model G: 
Vhot = 240 km s^l. Model J: Vhot = 100 km s'^. 



value as that in the case a^ot = 2 at the Milky Way scale, 
Vc = 220 kms-\ 

The CMRs in these models are shown in FigJlll The 
results for the model H with low feedback efficiency (right 
panels of Fig.[ll|) are almost the same as those for the model 
F (right panels of FigJq). So ahot does not affect the CMR 
when the feedback efficiency is low. On the other hand, in 
the model G with high feedback efficiency, the CMR is dif- 
ferent from that in the model E (left panels of Figs.FOl and 
H respectively). The slope of the CMR is flat at My ^ -21 
and steeper than that of observations at My ^ —21. 

In the case ahot ~ 5.5, the feedback strength /3 changes 
drastically at Vc ~ Vhot- So galaxies are divided into two 
classes at My ~ — 21. Probably this magnitude is given by 
Vhot = 240 km s~^ because it is almost the same circular 
velocity as that of the Milky Way galaxy, Vc — 220 km s~^, 
with My ~ —20.5. We find that the metallicity of the lu- 
minous class with Aly ^ —21 is nearly equal to 0.03, while 
that of the faint class with My <; —21 is nearly equal to 
0.015. Moreover, because the star formation in dwarf-scale 
galaxies is strongly suppressed, the age of galaxies becomes 
younger than that in the case of ahot ~ 2. Therefore the 
colours of dwarf ellipticals are bluer than the observations, 
while those of giant ellipticals are almost the same colour as 
the observations. Thus the slope of the CMR in the model 
G seems to be steeper than that of the observations. 



4.3 Effect of the UV background radiation 

In this subsection, we evaluate the effect of the Ultraviolet 
(UV) background radiation on the CMR. Now we obtain 
some evidences that the UV background radiation has ex- 



isted via the Gunn-Peterson test, the proximity effect, etc. 
The effect of the UV background on galaxy formation has 
been evaluated by Nagashima, Gouda & Sugiura (1999). 

The UV photons penetrate galactic gas clouds and ion- 
izing the envelopes of the clouds. We assume that absorption 
of the UV photons begins at the cooling radius rcooi and the 
UV photons are perfectly absorbed at a radius ruv- It is 
assumed that the hot gas which cannot cool does not affect 
the UV photons. The radius ruv is obtained by the inverse 
Stromgren sphere approximation, 



np{r)ne{r)a {Teq)4:nr dr 
= 7r(47rrcoo;)l-5 x 10''J-2i, 



(18) 



where Up and Ue are the proton and electron number den- 
sities respectively, a^^'{Teq) is the recombination coefficient 
to all excited levels at a temperature Teq, and J-21 is the 
normalized intensity of the UV radiation, J = J-21 x lO"'^^ 
erg cm~'^ s~^ sr~^ Hz~^. We solve this equation assuming 
the isothermal sphere distribution of gas in the clouds. The 
amount of newly cooled gas from hot gas is estimated as the 
amount of hot gas within ruv- 

We have many uncertainties about the evolution of the 
UV background, but at low redshift, it is suggested J_2i oc 
(1 + z)^ and J-2i{z = 2) ~ 1 (Pel 1995). At high redshift, 
there are too many uncertainties. So we approximate the 
UV background as J-2i{z = 2) = 1, J_2i oc (1 -I- 2)^,7 = 4 
at 2: < 2, 7 = — 1 at 2 < 2; < 5, and J-21 — at z > 5- 
Nagashima, Gouda & Sugiura (1999) showed that it is rea- 
sonable to adopt this type evolution of the UV background 
by comparing the luminosity function and the colour distri- 
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bution of galaxies given by their semi-analytic model with 
those given by observation. 

This effect on luminosity function of galaxies is similar 
to increasing feedback efficiency (see Cole et al. 1994; Na- 
gashima et al. 1999). Therefore, we introduce the effect of 
the UV background in our model instead of strengthening 
the feedback efficiency. 

In FigJl2, we show the CMRs in the models I and J, with 
short star formation time-scale, r ° = 2 Gyr. The feedback 
strengths are Vhot ~ 400 and the effect of the UV back- 
ground is not introduced in the model I, and Vhot = 280 
km s"'^ and the effect of the UV background is included in 
the model J. The observed slopes of the CMR are repro- 
duced well as well as the previous high feedback models. 

In both models, the slopes of the CMRs remain in the 
observational slope, while the colours becomes bluer. These 
show that the epoch of star formation delays by the strong 
feedback or the UV background because the star formation 
occurs in more massive galaxies by the effects. In these cases, 
we also confirmed that the metallicity-luminosity relations 
shows almost the same trend as the previous models. 

The virtue of the model J is that the amount of cold gas 
in 'Milky Way'-like galaxies, which belongs to a halo with 
circular velocity of 220 km s~^, is ~ lO^^M©. This value is 
nearly equal to that in our Galaxy. On the other hand, the 
model I {Vhot ~ 400 km s~^ and r ° = 2 Gyr), the amount 
of the cold gas ~ 2 x lO^M© is less than that in our Galaxy. 

Of course, such a degeneracy will be solved by observing 
other quantities, such as luminosity function, cold gas con- 
tent, and so on, because the dependences of these quantities 
on the feedback process, the star formation time-scale, and 
the UV background are different. This will be done in future 
work. 



5 CONCLUSIONS AND DISCUSSION 

We investigate the physical mechanisms determining the 
properties of the CMRs. We show that the feedback effi- 
ciency, the star formation time-scale, and the yield are the 
key parameters to determine the properties of the CMRs. 
Moreover, we find that the feedback process in the star- 
burst strongly affects the slope of the CMR when the star 
formation time-scale is long. 

The metallicity-luminosity relations, which are directly 
reflected to the CMRs, are formed as follows. At high lu- 
minosity, massive galaxies are dominated, whose feedback 
strength is weak. Because the feedback strength decreases as 
growing the mass of galaxies, the number of high metallicity 
stars increases. In this case, the mean stellar metallicity is 
determined by the yield mainly. 

On the other hand, at low luminosity, we must con- 
sider some mechanisms for suppressing the stellar metallic- 
ity. The following three parameters determine the mean stel- 
lar metallicity of the low luminosity ellipticals: the star for- 
mation time-scale, the feedback strength, and the intensity 
of the UV background [see eq.([lq)]. The final mean stellar 
metallicity increases when the star formation time-scale de- 
creases, because the time-scale of metal enrichment is given 
by the star formation time-scale and because the mean stel- 
lar metallicity increases as time passes monotonously in the 
case that a monolithic cloud forms stars. When the feed- 
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Figure 12. Same as FigM but for model I and J. The effect 
of the UV background radiation is considered in model J. Both 
models have a short star formation time-scale, T^ = 2 Gyr. Model 
I: Vhot = 400 km s^^ without the effect of the UV background. 
Model J: Vhot = 280 km s"! with the effect of the UV back- 
ground. 



back and UV background strength increase, the final mean 
stellar metallicity decreases. The reason for this is that the 
star formation is suppressed and then the amount of metals 
ejected into the gas decreases. 

Therefore one of the solutions to form the metallicity- 
luminosity relation is to increase the feedback strength when 
the star formation time-scale decreases. Another solution is 
to consider the effect of the UV background radiation. We 
find that the effect of increasing the feedback efficiency is the 
same effect as increasing the intensity of the UV background. 

The above conclusion suggests that the main reason for 
reproducing the CMR in the hierarchical clustering scenario 
is s imilar with that in the wind/collapse model (see Section 
3.3). That is, the feedback process of star formation is es- 
sential for making the metallicity-sequence and then, as a 
result, reproducing the slope of the CMR. Then from this 
point of view, we can say that the collapse/wind model suc- 
ceeds in showing the essential process for reproducing the 
CMR. However, in the traditional models such as the col- 
lapse/wind model, the formation epoch of elliptical galaxies 
is given by hand. In real, we should consider the cluster- 
ing process of dark haloes, cooling of baryonic gas, merging 
process of galaxies, etc. In fact, the age of elliptical galax- 
ies is not unique. In the semi-analytic model, we can pursue 
such formation processes of galaxies in the standard scenario 
based on the hierarchical clustering which is strongly sup- 
ported by the present observations of the large-scale struc- 
ture of the Universe. Moreover, in the hierarchical cluster- 
ing scenario, processes of suppressing star formation are not 
only the supernova feedback but also the UV background 
radiation. The former depends on the star formation time- 
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scale, while the latter affects all galaxies coeval. On the other 
hand, only the wind is the process of suppressing star forma- 
tion in the collapse/wind model. In this paper, we can show 
how such effects of suppression of star formation determine 
the slope of the CMR in the standard cosmological galaxy 
formation scenario. 

We also find that the feedback process in the starburst 
is important to make the CMR when the star formation 
time-scale is long enough, ~ 20 Gyr. When all of the cold 
gas and the metals do not escape from a galaxy but used to 
form stars in the starburst, this process is the same as the 
closed-box chemical enrichment model. If the star formation 
time-scale is long, the amount of the cold gas is large when 
the starburst occurs. Therefore the fraction of metal rich 
stars becomes large even in dwarf ellipticals. When the star 
formation time-scale is short enough (~ 2 Gyr), the feedback 
model in the starburst does not affect the CMR because 
almost all of the stars are formed in disc where the feedback 
is always effective. Such a value of star formation time-scale 
is often used for describing star formation in early-type disc 
galaxies (e.g., Arimoto, Yoshii & Takahara 1991). 

In our model based on the hierarchical clustering sce- 
nario, the metals released to hot gas return to cold gas later. 
We must evaluate whether this metal recycling affects the 
luminosity-metallicity relation. When t ^ t,, the reheated 
gas mass and newly formed stellar mass fractions are 



M, + Mreheat 
M, + Mreheat 



P + Rf 



l + /3-i?(l-/)^ 

1-R 
l + (3-R{l-f)- 



(19) 
(20) 



When considering small galaxies, which have a strong feed- 
back efficiency (/3 3> 1), most of cold gas are transformed to 
hot gas. So it is difficult for such small galaxies to form high 
metallicity stars. It means that the amount of metals formed 
in the galaxies is small. Therefore the hot gas does not evolve 
chemically. Thus the effect of metal recycling from hot gas 
to cold gas is slight. In order to confirm this, we calculate 
the CMR by using the model in which the metals returned 
from hot gas to cold gas are removed by hand. We find that 
the difference between the results of considering such met- 
als and removing metals is small, about 0.2 mag lower in 
V — K in all range of magnitude, and that the difference 
between the slopes is negligible. Therefore the effect of the 
metal recycling is negligible with regard to the slope of the 
CMR. 

In this paper, we show the physical mechanisms forming 
the CMRs and discuss the possible effects on the CMRs. 
We find a kind of 'degeneracy' among the star formation, 
the feedback, and the UV background for making the slope 
of the CMRs. This degeneracy will be solved by comparing 
other statistical properties of galaxies such as luminosity 
function and colour distribution in future work. Moreover, 
we should consider the tightness of CMRs. The tightness 
may reflect the dispersion of the mean stellar ages. However 
the colour of each galaxy is not simply estimated through its 
mean stellar metallicity and age. This is probably because a 
recent star formation by merging with small galaxies affects 
the colour. This effect is more emphasized in U-V rather 
than V-K, and the dispersions of the U-V CMRs are larger 
than those of the V-K CMRs especially in the case of r° — 20 



Gyr. We will investigate the origin of the tightness minutely 
in future work. 
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